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Abstract 

We set a preliminary 95% C.L. exclusion on the oscillation frequency of B° s — mixing 
using a sample of 400,000 hadronic Z° decays collected by the SLD experiment at the 
SLC between 1996 and 1998. The analyses determine the 6-hadron flavor at production 
by exploiting the large forward-backward asymmetry of polarized Z° — > bb decays as well 
as information from the hemisphere opposite that of the reconstructed B decay. In one 
analysis, B decay vertices are reconstructed inclusively with a topological technique, and 
separation between B® and B® decays exploits the B° s — > D~ cascade charge structure. In 
the other analysis, semileptonic decays are selected and the B decay point is reconstructed 
by intersecting a lepton track with the trajectory of a topologically reconstructed D 
meson. The two analyses are combined with a third analysis described elsewhere to 
exclude the following values of the B° s — B® mixing oscillation frequency: Am s < 7.6 ps _1 
and 11.8 < Am s < 14.8 ps _1 at the 95% confidence level. 
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1 Introduction 



Transitions between B° and B° mesons take place via second order weak interactions. In 
the Standard Model, a measurement of the oscillation frequency Arrid for B®-B® mix- 
ing determines, in principle, the value of the Cabibbo-Kobayashi-Maskawa matrix ele- 
ment \Vtd\ and constrains the Wolfenstein parameters p and (the CP- violating phase) rj, 
which are currently poorly constrained. However, theoretical uncertainties in calculating 
hadronic matrix elements are large (~ 20% Q]) and thus limit the current usefulness 
of precise Am^ measurements. Some of these uncertainties cancel when one considers 
the ratio between Am^ and Am s , leading to a reduced theoretical uncertainty (~ 4- 
8% m, U). Thus, combining measurements of the oscillation frequency of both B®-B® 
and B° s -B® mixing translates into a measurement of the ratio | V^I/jV^I and provides a 
stronger constraint on the parameters p and r\. 

Experimentally, a measurement of the time dependence of B°-B° mixing requires 
three ingredients: (i) the B decay proper time has to be reconstructed, (ii) the B flavor 
at production (initial state t — 0) needs to be determined, as well as (hi) the B flavor 
at decay (final state t = tdecay)- At SLD, the time dependence of B° s -B® mixing has 
been studied using three different methods, two of which are described below. The third 
method ("D s +tracks") is described elsewhere (see Ref. ||). All methods use the same 
initial state flavor tag but they use different techniques to reconstruct the B decay and 
tag its hnal state flavor. The data consists of some 400,000 hadronic Z° decays collected 
with the upgraded vertex detector (VXD3) between 1996 and 1998. The analyses exploit 
the large longitudinal polarization of the electron beam, P e = (73.4 ± 0.4)% for 1996-98, 
to enhance the initial state tag. 



2 Detector, Simulation and Event Selection 

The components of the SLD detector relevant to this analysis are presented here. The 
Liquid Argon Calorimeter (LAC) was used for triggering, event shape measurement and 
electron identification. It provides excellent solid-angle coverage (|cos#| < 0.84 and 
0.82 < | cos#| < 0.98 in the barrel and endcap regions, respectively). The LAC is divided 
longitudinally into electromagnetic and hadronic sections. The energy resolution for elec- 
tromagnetic showers is measured to be a / E = 15%/ ^E(GeV), whereas that for hadronic 

showers is estimated to be 60%/ J 1 E{GeV). The Warm Iron Calorimeter (WIC) provides 
efficient muon identification for | cos#| < 0.60. Tracking is provided by the Central Drift 
Chamber (CDC)@] for charged track reconstruction and momentum measurement and 
the CCD pixel Vertex Detector (VXD)|5|] for precise position measurements near the in- 
teraction point. Aside from the WIC, these systems are immersed in the 0.6 T field of the 
SLD solenoid. Charged tracks reconstructed in the CDC are linked with pixel clusters in 
the VXD by extrapolating each track and selecting the best set of associated clusters||]. 
The track impact parameter resolutions at high momenta are 7.8 pm and 9.7 pm in the 
r0 and rz projections respectively (z points along the beam direction), while multiple 
scattering contributions are 33 pm /(psm 3 ^ 2 6) in both projections (where the momentum 
p is expressed in GeV/c). 
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The centroid of the micron-sized SLC Interaction Point (IP) in the rip plane is 
reconstructed with a measured precision of ajp = (4 ± 2) /zm using tracks in sets of ~ 30 
sequential hadronic Z° decays. The median z position of tracks at their point of closest 
approach to the IP in the r<p plane is used to determine the z position of the Z° primary 
vertex on an event-by-event basis. A precision of ~ 20 fim on this quantity is estimated 
using the Z° — > bb Monte Carlo (MC) simulation. 

The simulated Z° — > qq events are generated using JETSET 7.4 ||. The B meson 
decays are simulated using the CLEO B decay model [7| tuned to reproduce the spectra 
and multiplicities of charmed hadrons, pions, kaons, protons and leptons as measured at 
the T(4S) by ARGUS and CLEO §. Semileptonic decays of B mesons follow the ISGW 
model including 23% D** production. The branching fractions of the charmed hadrons 
are tuned to the existing measurements [[HJ. The lifetimes of B mesons and 6-baryons 
are chosen to be t b + = 1.656 ps, t b q = 1.562 ps, t b o = 1.464 ps, and T\ b = 1.208 ps, 
according to recent world averages [pT| . The 6-quark fragmentation follows the Peterson 



et al. parameterization [E|. Finally, the SLD detector is simulated using GEANT 3.21 
0- 

Hadronic Z° event selection requires at least 7 CDC tracks which pass within 5 cm 
of the IP in z at the point of closest approach to the beam and which have momentum 
transverse to the beam direction p±_ >200 MeV/c. The sum of the energy of the charged 
tracks passing these cuts must be greater than 18 GeV. These requirements remove back- 
ground from Z° — > events and two-photon interactions. In addition, the thrust axis 
determined from energy clusters in the calorimeter must have |cos#r| < 0.85, within 
the acceptance of the vertex detector. These requirements yield a sample of ~ 311,000 
hadronic Z° decays. 

Good quality tracks used for vertex finding must have at least two associated VXD 
hits and p± >250 MeV/c. Additional requirements are that the tracks have either three or 
more VXD hits or else satisfy the following criteria: (i) have a CDC hit at a radius<50 cm, 
(ii) have >23 hits to insure that the lever arm provided by the CDC is appreciable, (iii) 
extrapolate to within 1 cm of the IP in rcf) and within 1.5 cm in z to eliminate tracks 
which arise from interaction with the detector material, (iv) have a x 2 /d.o.f.< 8 for both 
the CDC portion of the track and the combined VXD-CDC track. 

Tracks reconstructed in the vertex detector but unsuccessfully tracked in the drift 
chamber are also used in the analyses. Such "VXD-only" tracks are constructed from hits 
in all three of the vertex detector layers and are used primarily to improve the overall 
vertex charge reconstruction. 

Both analyses make use of the inclusive topological vertexing technique |14[] devel- 
oped for B lifetime |T5j and Rb [0] analyses to tag and reconstruct 6-hadron decays. The 
b purity of the sample is increased by reconstructing the vertex mass M, which includes 
a partial correction for missing decay products (see Ref. This inclusive vertexing 

technique has been adapted for semileptonic decays to reconstruct the D decay topology 
(see below). 
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3 Initial State Flavor Tagging 



The large forward-backward asymmetry for polarized Z° — > bb decays is used as a tag of 
the initial state flavor. The polarized forward-backward asymmetry A FB can be described 
by 

Afb (cos 9t) = 2A b — — — — , 1 

1 — A e P e 1 + cos 2 9t 

where A b = 0.935 and A e = 0.150 (Standard Model values), P e is the electron beam 
longitudinal polarization, and &t is the angle between the thrust axis and the electron 
beam direction (the thrust axis is signed such that it points in the same hemisphere as 
the reconstructed B vertex). Thus, left- (right-) polarized electrons tag b (b) quarks in the 
forward hemisphere, and b (b) quarks in the backward hemisphere. Averaged over our 
acceptance, this yields an average correct tag probability of 0.74 for an average electron 
polarization P e = 73%. The probability for correctly tagging a b quark at production is 
expressed as 

P A (cos6 T ) = . (2) 

A jet charge technique is used in addition to the polarized forward-backward asym- 
metry. For this tag, tracks in the hemisphere opposite that of the reconstructed vertex 
are selected. These tracks are required to have momentum transverse to the beam axis 
p± > 0.15 GeV/c, total momentum p < 50 GeV/c, impact parameter in the plane per- 
pendicular to the beam axis 5 < 2 cm, distance between the primary vertex and the track 
at the point of closest approach along the beam axis Az < 10 cm, and | cos6*| < 0.90. 
With these tracks, an opposite hemisphere momentum-weighted track charge is defined 
as 

Qopp = ^qi \pi ■ t\ , (3) 

i 

where is the electric charge of track i, f>i its momentum vector, T is the thrust axis 
direction, and k is a coefficient chosen to be 0.5 to maximize the separation between b 
and b quarks. The probability for correctly tagging a b quark in the initial state of the 
vertex hemisphere can be parameterized as 

P ^Qo PP ) = TT ^ ; , (4) 

where the coefficient a = —0.27, as determined using the Monte Carlo simulation. This 
technique yields an average correct tag probability of 0.66 and is independent of the 
polarized forward-backward asymmetry tag. 

Finally, the tag is further enhanced by the addition of other flavor-sensitive quan- 
tities from the hemisphere opposite that of the selected vertex. For this purpose, the in- 
clusive topological vertexing technique mentioned earlier is used. The sensitive variables 
are: the total track charge and charge dipole of a topologically reconstructed vertex, the 
charge of a kaon identified in the Cherenkov Ring Imaging Detector, and the charge of a 
lepton with high transverse momentum with respect to the direction of the nearest jet. 
The addition of these tags improves the average correct tag probability by about 0.03. 
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The various tags are combined, taking correlations into account, to form an over- 
all initial state tag characterized by a 6-quark probability Pj. The average correct tag 
probability is 0.75 for the lepton+D and 0.77 for the charge dipole analyses with 100% 
efficiency Fig. [I] shows the Pi distributions for data and Monte Carlo in the charge 
dipole analysis (see below for a description of the analysis), and also indicates the clear 
separation between b and b quarks. 




ii= i i i I i i i I i i i I i i i I i i i ^~ 
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Initial State b-auark Probability 

Figure 1: Distribution of the computed initial state b-quark probability for data (points) 
and Monte Carlo (histograms ) showing the b and b components for the events selected in 
the charge dipole analysis. 



4 Lepton+D Analysis 

The lepton+ "D" analysis aims at reconstructing the B and D vertex topologies of semilep- 
tonic B decays. It proceeds by first selecting event hemispheres containing an identified 
lepton (e or jj) with | cos^l < 0.7. Then, a D vertex candidate is reconstructed using the 
inclusive topological technique described earlier. If multiple displaced vertices are found 
in the same hemisphere, the vertex with the largest invariant mass is chosen to be the D 
vertex candidate. A resultant D "track" is created using the D vertex location and the 
parameters of all tracks attached to it. Furthermore, the D track error matrix is corrected 
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0.2 0.4 0.6 0.8 1 

Neural Network Output 

Figure 2: Distribution of neural network output variable for data (solid points) and Monte 
Carlo (solid histogram). Also shown are the Monte Carlo distributions for B® right-sign 
leptons (open circles), B% right-sign leptons (dashed histogram) and B® wrong-sign leptons 
(dotted histogram). 

to take into account the fact that the D decay is not fully reconstructed. Finally, the B 
decay vertex is reconstructed by intersecting the lepton and D tracks. 

A neural network is used to clean up the D vertex candidates and reduce the 
contamination from cascade (b — > c — > I) charm semileptonic decays. The Jetnet neural 
network package is used with 6 inputs and 12 hidden nodes. The inputs are the transverse 
momentum of the lepton with respect to the B vertex direction (vector stretching from 
the IP to the B vertex), the B decay length (magnitude of that vector), the transverse 
momentum of the lepton with respect to the D vertex direction (vector stretching from 
the B vertex to the D vertex), the mass M of the charged tracks associated with the B 
decay and the distance of closest approach of the lepton to the B vertex. A distribution 
of the neural network output is shown in Fig. ||]. A minimum cut of 0.65 is applied on the 
neural network output. 

For this analysis, only vertices with positive reconstructed decay length are se- 
lected. To enhance the fraction of B° s decays, the sum of lepton + D vertex track charges 
is required to be Q = 0. This enhances the B° s fraction to 16.4% of all b hadrons in 
the Z° bb MC (the B° s production fraction in the Z° -> bb MC is 10.0%). The udsc 
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Figure 3: Distributions of lepton momentum transverse to the D vertex trajectory, D 
vertex track multiplicity, D vertex mass and lepton+D vertex mass for data (points) and 
Monte Carlo (histograms). 



contamination is 0.5% in the final sample. 

The B° s fraction is further enhanced by using the Cherenkov Ring Imaging Detector 
to identify kaon candidates among the D vertex tracks. For the subsample containing a 
lepton-kaon pair with opposite charge, the fraction increases to 38.5% and the udsc 
contamination remains at 0.5%, as determined from MC. 

A sample of 2087 decays is thus obtained in the 1996-98 data. Various comparisons 
between data and Monte Carlo simulation were performed which showed good agreement. 
For example, Fig. ^ shows the distributions of lepton momentum transverse to the D 
vertex trajectory, D vertex track multiplicity, and invariant mass of all tracks in the D 
vertex and (assuming all tracks are pions) as well as in both B and D vertices. 

A powerful check of the analysis and the purity of the final state tag is the 
polarization-dependent forward-backward asymmetry shown in Fig. f|. A clear asym- 
metry is observed, in reasonable agreement with the Monte Carlo, indicating that the 
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SLD Lepton+D 
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Figure 4: Distribution ofcosO for the thrust axis direction signed by the product (Qi ep tc 
P e ) for data (points) and Monte Carlo (histogram) . 
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Figure 5: Distributions of the decay length and relative boost residuals for B® 
decays in the simulation. 

8 



SLD lepton+D 
900 




Figure 6: Distributions of the fraction of decays tagged as "mixed" as a function of 
reconstructed proper time (left) and reconstructed proper time (right) for the data (points) 
and the likelihood function (histograms). 



final state tag purity is adequately modeled in the simulation. 

The study of the time dependence of B®-B® mixing requires a precise determination 
of the B decay proper time t = L/(7/3c), where L is the reconstructed decay length 
(distance between the IP and the B vertex) and 7/3 = Pb/^ib is computed from the 
estimated B momentum ps and the known mass of the B meson, m#. Reconstruction 
of the 6-hadron boost uses both tracking and calorimeter information. A description of 
the reconstruction algorithm may be found in Ref. jfT/j] . The overall performance of the 
decay length and boost measurements for B° s decays proceeding via the direct (b — > /) 
transition is shown in Fig. [5| The proper time distribution is shown in Fig. ^. 

The final state B° or B° flavor is tagged by the sign of the lepton charge. Each 
decay is assigned a final state 6-quark probability P/, defined such that Pf > 0.5 (< 0.5) 
corresponds to a negatively (positively) charged lepton which then tags the decay as B 
(B). The magnitude of the correct tag probability depends on the sample composition 
as well as on the lepton p?. The lepton sources in selected B° s decays are as follows: 
94.8% (b -> /-), 2.6% (b -> c -> l+), 1.4% (b -> c -> 0.6% (b -> X~) (right-sign 
misidentified lepton), and 0.6% (6 — > X + ) (wrong-sign misidentified lepton). The final 
state correct tag probability is thus 0.968. Further enhancement of the tag is achieved by 
taking into account the strong p? dependence of the various lepton source fractions. 
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4.1 Likelihood Function 



The search for the time dependence of B° s -B® mixing is carried out with a likelihood 
analysis which includes the effect of detector smearing, mistag of both initial and final 
states, selection efficiencies and the dependence on the oscillation frequency Am s . The 
probability that a meson created as a B® (B®) will decay as a B® (B®) after proper time 
t can be written as 

P u (t) = ^e- rt [1 + cos(Am s t)} , (5) 

where Am s is the mass difference between the mass eigenstates, T is the average decay 
width of the two states and P u denotes the probability to remain 'unmixed'. The effects 
of CP violation are assumed to be small and are neglected. Similarly, the probability that 
the same initial state will 'mix' and decay as its antiparticle is 

P m (t) = \z~ Tt [l-cos(Am,f)] . (6) 

Decays are tagged as mixed or unmixed if the product (Pi — 0.5) x (Pf — 0.5) is 
smaller or greater than 0, respectively. The probability for a decay to be in the mixed 
sample is expressed as: 

e -t/r u 



V mixed (t, Am s ) =/„ Wbfi + g c ; + g x J) + (1 - + g x u w ) 

(9i l + 97 + + [2^ 7 - 1] cos Am d t) + (g™ + g™)(l - [2w< - 1] cos Am d t)_ 



+ fde 



-t/r d 



2 r d 

t/r s 



+ f — 1 \(gf + g c I + + [2w r - 1] cos Am s t) + (g™ + g™)(l - [2W 1 - 1] cos Am s t) 

2 r, L 



is 

+ fbaryon—- W (gbaryon + dbaryon + dbaryon) + — W )(9baryon + 9baryon) 

'baryon 

+ ^ F udsc (t), (7) 

where fj represents the fraction of each 6-hadron type and background (j = u, d, s, A, 
and udsc correspond to B + , B®, B®, 6-baryon, and udsc background), Tj is the lifetime 
for b hadrons of type j, w 1 is the initial state mistag probability, gf, gj and g x - r are the 
fractions of right-sign (b — > l~), (b — > c — > l~) and (6 — > X~) leptons, respectively, g c ™ 
and (yfj 1 " are the fractions of wrong-sign (b — > c — > / + ) and (6 — > X + ) leptons, respectively, 
and F udsc (t) is a function describing the proper time distribution of the udsc background 
(a sum of two exponentials is used). A similar expression for the probability V unm i xed to 
observe a decay tagged as unmixed is obtained by replacing the mistag rate w 1 by 1 — w 1 . 

Detector and vertex selection effects are introduced by convoluting the above prob- 
ability functions with a proper time resolution function 1Z(T, t) and a time-dependent 
efficiency function e(t): 

P mixed (T, Am s ) = / V mixed (t, Am s ) K(T, t) e(t) dt , (8) 
Jo 
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where t is the "true" time and T is the reconstructed time. Again, a similar expression 
applies to the unmixed probability Punmixed- The resolution function is parameterized by 
the sum of four Gaussians: 



*{T,t) = E E k 7Tw^ e H=W , (9) 
»=i i=i <Tij[t)y/ 2n 

where the index i — 1 (i = 2) corresponds to the core (tail) component of the decay 
length resolution an, the index j = 1 (j = 2) similarly corresponds to the core (tail) 
component of the relative boost resolution a^pj/^p. The various fractions are fn = 0.36, 
fn = /21 = 0.24, and f 2 2 = 0.16. The proper time resolution <Jij(t) is a function of 
proper time that depends on the measured boost 7/?, its resolution and the decay length 
resolution: 



<7«(t) 



1/2 



(10) 



For each decay, the resolution <jl is computed from the vertex fit and IP position mea- 
surement errors, with a scale factor determined using the MC simulation (the scale factor 
is introduced mostly to account for the fact that the analysis does not attempt to fully 
reconstruct the D meson decay). The relative boost residual o^l^fi is parameterized as a 
function of the lepton + D vertex total track energy, with parameters extracted from the 
MC simulation. Similarly, the efficiency e(t) is parameterized using the MC simulation. 
All parameterizations are performed separately for each 6-hadron type. For example, the 
efficiency for B° s decays is given by 

1 - e bt 

with a = 0.148, b = —5.7, and c = 0.0072. Furthermore, and a 7j g resolutions are 
handled separately for the main lepton sources (&—►/), (b — > c(c) — ► /) and (b — > X). 
As a consequence, different resolution functions are used for the different sources and the 
expressions for P m i xe d and P unm ixed are modified accordingly. 

The likelihood function is constructed from the calculated probabilities for events 
tagged as mixed and unmixed. The total likelihood for the sample is given by 

#mixed ^unmixed 

£ W Pmixed J^J Punmixed ■ (12) 

i=l j=l 



4.2 Oscillation Analysis 

The study of the time dependence of B®-B® mixing is carried out using the amplitude 
method described in Ref. ||18|| . Instead of fitting for Am s directly, the analysis is performed 
at fixed values of Am s and a fit to the amplitude A of the oscillation is performed, i.e. 
in the expression for the mixed and unmixed probabilities, one replaces [1 ± cos(Am s t)} 
with [1 ± Acos(Am s t)]. This method is similar to Fourier transform analysis and has the 
advantage of facilitating the combination of results from different analyses and different 
experiments. 
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The measured amplitude for the lepton+D analysis is shown as a function of Am s 
in Fig. [7[ The measured values are consistent with A = for the whole range of Am s up 
to 30 ps -1 and no evidence is found for a preferred mixing frequency. 

Systematic uncertainties have been computed following Ref. jTj| and are summa- 
rized in Table |I] for several Am s values. Uncertainties in the sample composition are 
estimated by varying the fraction of udsc background by ±20% and the production frac- 
tions of B® and 6-baryons according to 0.100 ± 0.012 and 0.099 ± 0.017, respectively JTT|. 
Other physics modeling uncertainties are t(B + ) = 1.656±0.025 ps, t(B®) = 1.562±0.029 
ps, r(B° s ) = 1.464 ± 0.057 ps, r(A & ) = 1.208 ± 0.051 ps, and Am d = 0.480 ± 0.020 ps -1 . 
Uncertainties in the modeling of the detector include ±10% and ±20% variations in decay 
length and boost resolutions, respectively. Initial state tag uncertainties are estimated by 
varying the correct tag probability by ±0.02 (i.e., a ±10% variation of the mistag rate), 
corresponding to the expected contribution from uncertainties in the measured electron 
beam polarization, the value of Ab, and the self-calibrated jet charge analyzing power. Fi- 
nal state tag uncertainties include a ±15% variation in the lepton misidentification rate, 
as well as the effect of uncertainties in the branching ratios B(b — > I) = 0.112 ± 0.002, 
B{b -> c -> I) = 0.016 ± 0.004, and B(b -»• c Z) = 0.080 ± 0.004. The dominant 
uncertainty is the B® production fraction in Z Q — ► bb events. 
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Table 1: Measured values of the oscillation amplitude A with a breakdown of the main 
systematic uncertainties for several Am s values in the lepton+D analysis. 

Am s 

Measured amplitude A 



a 



syst 



5 ps _1 


10 ps" 1 


15 ps _1 


20 ps~ 


-0.644 


0.433 


0.668 


2.080 


±0.482 


±0.714 


±1.070 


±1.310 


+0.113 


+0.503 


+0.249 


+0.292 


g 442 


n 253 


—0 389 


n ]fig 


— U.134 


— 0.U76 


— U.lby 


— U.U2V 




in onfi 
-ru.zyo 


in 991 


in 9fl9 


+0.047 


+0.065 


+0.042 


+0.025 


-0.051 


-0.096 


-0.041 


-0.027 


+0.003 


-0.007 


+0.028 


+0.005 


-0.002 


+0.008 


-0.033 


-0.044 


-0.162 


+0.362 


+0.012 


+0.206 


-0.040 


-0.113 


-0.327 


-0.138 


+0.033 


+0.049 


+0.033 


+0.002 


-0.061 


-0.176 


-0.034 


-0.036 


-0.093 


+0.010 


+0.000 


+0.004 


-0.008 


-0.028 


-0.000 


-0.005 


+0.014 


-0.001 


-0.022 


-0.025 


-0.226 


+0.003 


+0.029 


+0.017 


+0.012 


+0.163 


+0.064 


+0.019 


-0.195 


-0.000 


-0.070 


-0.024 


+0.008 


+0.015 


-0.005 


+0.000 


-0.010 


-0.025 


+0.005 


-0.004 



T. = B(b =; b°) 

/a = B(b — > 6— baryon) 
decay length resolution 
boost resolution 
B° s lifetime 
Am d 

initial state tag 

B(6->Z), B(6->c->Z), B{b 

lepton misidentification 







5 Vertex Charge Dipole Analysis 

The charge dipole analysis aims at selecting decays with distinct B and D vertices and 
tags the B° or B° decay flavor based on the charge difference between them. This analysis 
technique was first developed by SLD and relies extensively on the excellent resolution of 
the vertex detector. 

In the following, we first describe the algorithm used to identify primary, secondary 
and tertiary vertices, then discuss details of the B®-B® mixing analysis. 

5.1 Ghost Track Algorithm 

The B decay flavor tag with the charge dipole relies on the kinematic fact that the boost 
of the B decay system carries the cascade charm decay downstream from the B decay 
vertex. Monte Carlo studies show that in B decays producing a single D meson the 
cascade D decays on average 4200 /xm from the IP, while the intermediate B vertex is 
displaced on average only 46 transversely from the line joining the IP to the D decay 
vertex. This kinematic stretching of the B decay chain into an approximately straight 
line is exploited by the ghost track algorithm. This new algorithm has two stages and 
operates on a given set of selected tracks in a jet or hemisphere. First, the best estimate 
of the straight line from the IP directed along the B decay chain is found. This line is 
promoted to the status of a track by assigning it a finite width. This new track, regarded 
as the resurrected image of the deceased B hadron, is called the "ghost" track. Secondly, 
the selected tracks are vertexed with the ghost track and the IP to build up the decay 
chain along the ghost direction. Both stages are now described in more detail. 
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Given a set of tracks in a hadronic jet or hemisphere a new track G is created with 
the properties that it is a straight line from the IP directed along the jet or thrust axis 
and has a constant resolution width of 25 /mi in both r<p and rz. For each track i a vertex 
is formed with track G and the vertex location ri, fit x\ and Lj are determined (Lj is the 
longitudinal displacement from the IP of ri projected onto the direction of track G). This 
is calculated for each of the tracks and the summed x 2 is formed: 

y 2 - V x * U ~ °-° (13) 

X * ~ ^ (2x1 - X 2 ) U < 0.0 {l6) 

where Xoi is the xt °f track % to track G determined at Lj = rather than at the best fit 
vertex location. The aim is to construct this quantity, x\i such that when the direction 
of G is varied the minimum of x| provides the best estimate of the B decay direction. If 
the initial direction is a relatively long way from the B line of flight, some or all of the 
decay tracks may vertex with G with a negative value of Lj. In this case the 2%^ — x 2 
term above helps to push track G towards the B flight path as Xs * s minimized. This first 



minimization using equation [13| is designed for this purpose. (Note that the contribution 
of each track as x"s i s minimized changes in a continuous manner even if Lj changes sign 
since x\ = Xm at U = 0.) 

The value of x| is recalculated as track G is rotated (about the pivot at the IP) 
incrementally in ever decreasing angular steps 58 and 5<p until the minimum is found 
within the required precision (< 0.1 mrad, i.e. within 1 fim at 1 cm from the IP). The 
width of track G is set such that the maximum x\ = 1-0 f° r all tracks with Lj > (if 
this is less than 25/im, it is restored to 25/xm). The track G is now consistent with all 
potential B decay candidate tracks (Lj > 0) at the level Xi < 1-0- In other words, the 
new width of G measures the degree to which the tracks conform to a straight line decay 
chain. A second iteration in 58,5<j) now takes place with the summed x 2 redefined as: 

which is not sensitive to any spurious background track with a negative value of Lj which 
might otherwise perturb the direction of track G. After finding the new minimum of Xs 
the width of G is again recalculated such that < 1.0 for all tracks i with Lj > 0. Again 
this width is required to be at least 25 /xm. Track G is now directed along the best guess 
of the B decay line of flight and has a width such that it is consistent with potential B 
decay tracks in the jet, track G is now called the "ghost" track. 

The second stage of the algorithm begins by defining a fit probability for a set of 
tracks to form a vertex with each other and with the ghost track (or IP). This probability 
then measures the likelihood of the set of tracks both belonging to a common vertex and 
being consistent with the ghost track (or IP) and hence forming a part of the B decay 
chain. These probabilities are determined from the fit x 2 which is in turn determined 
algebraically from the parameters of the selected tracks and the ghost track (or the 7 x 
7 x 30/im 3 ellipsoid assumed for the IP). The earlier requirement that each Lj > track 
makes a x\ ^ 1-0 with the ghost track has the effect that the fit probabilities have the 
desired property of having an approximately flat distribution from 0.0 to 1.0 for genuine 
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vertices, independent of both multiplicity and decay length. This property also relies on 
the choice of the number of degrees of freedom as 2N— 2 (or 2N) when fitting N tracks 
together with the ghost track (or IP). Fake vertices peak at probability close to 0.0. 

For a set of N tracks, there are initially N+l candidate vertices (N 1-prong sec- 
ondary vertices and a bare IP). A matrix of track % - track j associations is constructed 
to store the calculated probabilities of each candidate vertex pair fitted together with the 
ghost track. A further column and row is added to the matrix to store the probabilities of 
each track fit with the IP ellipsoid. The upper triangle of the matrix (i.e. the ij (i < j) 
elements) stores the probabilities while the lower triangle (initialized with ij (i > j) ele- 
ments set to 0.0) indelibly records which tracks (and IP) have been assigned together in a 
common vertices as the algorithm progresses. Once the upper triangle has been filled, the 
highest probability in the matrix table is found and the corresponding candidate vertex 
pair are from then on tied together in a new candidate vertex for all future computations 
by flagging the corresponding lower triangle elements of the matrix with non-zero values. 
The upper triangle of the matrix is now refilled taking into account the associations that 
have so far been made, the new maximum probability is found, and the corresponding 
subset of the tracks and IP is tied together. At each iteration of combining the maximum 
probability matrix element contributors, the number of candidate vertices decreases by 
one. The iterations continue until the maximum probability is less than 1%. At this 
point the tracks and IP have been divided into unique subsets by the associations thereby 
defining topological vertices. 

Jets or hemispheres in which three vertices are found - the primary (which includes 
by definition the IP), a secondary and a tertiary - are used for the charge dipole analysis. 
The secondary vertex is identified as the B decay vertex and the tertiary as the cascade 
charm decay. As well as improving the purity and efficiency of the dipole reconstruc- 
tion (by requiring the vertices be consistent with a single line of flight) the ghost track 
algorithm has the additional advantage of allowing the direct reconstruction of 1-prong 
vertices, including the topology consisting of 1-prong B and D decays. 

Finally, VXD-only tracks are attached to the B decay chain to improve the overall 
charge reconstruction but they are not used in the determination of the B and D vertex 
positions. The attachement criteria rely on variables % and Lj for each track i, as defined 
below. A vertex axis is formed by a straight line joining the IP to a vertex combining 
both secondary and tertiary tracks. The 3-D distance of closest approach of the track 
to the vertex axis, T;, and the distance from the IP along the vertex axis to this point, 
Li, are calculated. Track % is attached to the vertex if % < 0.1 cm, Li > 0.025 cm and 
0.25 < Li/ L < 2 (where L is the distance between the IP and the combined vertex). An 
average of 0.2 VXD-only tracks is added per decay, to be compared with an average of 
5.0 VXD+CDC tracks per decay (for M > 2 GeV/c 2 ). The VXD-only tracks attached 
to the B decay chain are further attached to either the B or D vertex according to their 
longitudinal displacement Lf tracks with Li < Lb + 0.5(Ld — Lb) are attached to the B 
vertex and all others are attached to the D vertex, Lb {Ld) is the distance between the 
IP and the B (D) vertex. 
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5.2 Event Selection 



Hemispheres containing both a secondary and a tertiary vertex are selected for the charge 
dipole analysis. Furthermore, the invariant mass computed using all secondary and ter- 
tiary vertex tracks is required to be M > 2 GeV/c 2 (the computed mass includes a 
partial correction for missing decay products) and the total track charge Q to t (from both 
secondary and tertiary vertices) is required to be zero to enhance the fraction of B° s decays 
in the sample and to increase the quality of the charge difference reconstruction for neutral 
B decays. As mentioned in the previous section, the (secondary) vertex that is closer to 
the IP is labelled "JB" and that further away (tertiary) is labelled "D." A "charge dipole" 
is defined as SQ = Dbd x SIGN(Qd — Qb), where Dbd is the distance between the two 
vertices and Qb (Qd) is the charge of the B (D) vertex. Positive (negative) values of 
SQ tag B° (B°) decays. Requirements on the vertices are: 250 /im < D B d < 1 cm , D 
vertex mass < 2.0 GeV/c 2 (assuming all tracks are pions), B vertex decay length L > 0, 
Qb 7^ Qd, "ghost" track width < 300 /im and cosine between the straight line connecting 
the IP and the B vertex, and the nearest jet axis direction < 0.9. The decay is rejected 
if any attached VXD-only track has pj_ > 4 GeV/c, since in that case the charge is not 
reliably reconstructed. MC studies indicate that, after these selection cuts, the track as- 
signment to the B (D) vertex is 84% (86%) correct for B® decays containing one D meson 
in the final state, i.e. 84% (86%) of all tracks in the B (D) vertex originate from the decay 
point of the B (D) meson. For all data and MC events, hemispheres already containing a 
vertex selected by the lepton+D or D s +track analyses are removed to keep the analyses 
statistically uncorrelated. The udsc background is further suppressed by demanding that 
the event contains either an opposite hemisphere topological vertex with M > 2 GeV/c 2 
or at least 3 tracks with positive 2-D impact parameter > 3 o. The udsc fraction is thus 
reduced to 2.2%. 

Applying all the above cuts, a sample of 8556 decays is selected in the 1996-98 
data. Figure |8| shows distributions of the B and D vertex track multiplicities, as well as 
the distance and charge difference between B and D vertices in the selected sample. Good 
agreement between data and MC is obtained. A slight discrepancy in the D vertex track 
multiplicity is apparent but was determined to have negligible impact on the analysis. 
Requiring that the total track charge be zero boosts the B° s fraction from its assumed 
production value of 10.0% to 15.3%. Figure |9| displays the distribution of charge dipole 
SQ for the data sample and also indicates the separation between b hadrons containing b 
or b quarks in the MC. 

The average correct tag probability for the charge dipole tag is 0.76 for selected B° s 
decays and is parameterized as a function of decay length, as shown in Fig. [10]. Further- 
more, the correct tag probability depends on the charm content of the decay products. 
Thus, the correct tag probability is parameterized separately for B® and B® decays into 
five different final states: D°X, D + X, D S X, charmed hadron X, and DDX (this last 
category also incorporates charmonium production, i.e. it includes all b — > ccs decays). 
For example, the correct tag probability is 0.88 for B® — ► D S X decays but only 0.53 for 
B® — ► DDX decays. The above correct tag probabilities are extracted from the MC 
simulation and include a multiplicative dilution factor Sp = 0.95. This factor is conser- 
vatively assigned to reproduce the overall fraction of decays tagged as mixed observed in 
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Figure 8: Distributions of B and D vertex track multiplicity, as well as distance and 
charge difference between B and D vertices for data (points) and Monte Carlo (histograms) 
in the charge dipole analysis. 
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Figure 9: Distribution of the charge dipole for data (points) and Monte Carlo (solid 
histogram). Also shown are the contributions from b hadrons containing a b quark (dotted 
histogram) or a b quark (dashed histogram). 
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Figure 10: Charge dipole correct tag probability as a function of reconstructed decay 
length in simulated B° s — > D S X (solid circles) and B® — > DDX (open squares) decays. 
The functions are the fit results used to parametrize the charge dipole correct tag probability 
as a function of decay length. 
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the data, as well as to properly describe the time evolution of B®-B® mixing in the data 
(see below). 

As hadronic decays of B mesons are not as well known as semileptonic decays, it is 
important to check the correct tag probability estimated using measured quantities like the 
polarization-dependent forward-backward asymmetry shown in Fig. [11]. Good agreement 
with the MC is observed, indicating that the charge dipole correct tag probability is well 
modeled. It should be noted that this asymmetry is diluted by both initial and final state 
mistags and by B°-B° mixing. The dilution due to mixing can be reduced by selecting 
vertices with total charge Qtot = ±1, in which case a stronger asymmetry is observed (see 
Fig. |TTJ) . Another useful test of the charge dipole tag in B% decays is the measurement 
of the time dependence of B^-B^ mixing. This has been checked using the full likelihood 
analysis described in the following section. Fitting for the B®-B® mixing frequency yields 
Am<i = 0.495 ± 0.032 ps _1 (statistical error only), see Fig. [12|. This value agrees well 



with the latest world average value of 0.485 ± 0.015 ps 1 [fT9|. The mixed fraction as a 




function of proper time is displayed in Fig. [13]. The figure shows the expected increase in 
the fraction of decays tagged as mixed for decays at small proper time. Most of this effect 
originates from misreconstructed B + decays near the IP which tend to have random final 
state tags. 



5.3 Likelihood Function 

The B° s -B® mixing fit for the charge dipole analysis is performed in a way similar to the 
lepton+D analysis. Decays are tagged as mixed or unmixed if the product (Pj — 0.5) x 
(Pf — 0.5) is smaller or greater than 0, respectively. The probability for a decay to be in 
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where /j represents the fraction of each 6-hadron type and background (j = u, d, s, baryon, 
and udsc correspond to B + , B d , B®, 6-baryon, and udsc background), F udsc {t) is a function 
describing the proper time distribution of the udsc background (a sum of two exponentials 
is used), Tj is the lifetime for b hadrons of type j, w 1 is the initial state mistag probability, 



wf, and w^ aryon 



are the final state mistag probabilities for B + and b baryons, whereas w t 



and Wg k are the final state mistag probabilities for B d and B®, with the index k = 1, 5 
representing the five different decay final states: D°X, D + X, D S X, charmed hadron X, 
and DDX, g d k and g sk are the fractions of B d and B® decays into each of the above final 
states. A similar expression for the probability V un mixed to observe a decay tagged as 
unmixed is obtained by replacing the initial state mistag rate w 1 by (1 — w 1 ). 

Several of the quantities in Eq. (|15|) are determined on an event by event basis. The 
initial state mistag probability w 1 depends on cos 9 of the thrust axis, the electron beam 
polarization, as well as several quantities from the opposite hemisphere: jet charge, vertex 
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charge, kaon charge, lepton charge and dipole charge. The final state mistag probabilities 
Wp. depend on the reconstructed decay length to take into account the degradation of 
the charge dipole tag close to the IP. This effect is fairly weak for B%, B® and b baryon 



decays (see, for example, Fig. [T/y) but it is significant for B + decays. Finally, the decay 
final state fractions g^k and g s k are parametrized as a function of the overall vertex mass 
M. For example, the fraction of B® decays into DDX final states decreases from 0.28 at 
M = 2 GeV/c 2 to 0.12 at M = 5 GeV/c 2 . 

As described in Sec. |4.1| , the functions V m i xe d and V unm i xe d are convoluted with 
a proper time resolution function TZ(T,t), see Eq. (|) and a time- dependent efficiency 
function £j(t). Separate efficiency functions are extracted for each b hadron type using 
the simulation. The relative boost resolution <7 7j g / 7/3 is parametrized with the sum of two 
Gaussians using the MC simulation. Considering all selected B® decays, the widths of 
the two Gaussians are o m = 0.07 and <tb 2 = 0.21, where the first Gaussian represents 
60% of the decays. However, the analysis takes into account the strong dependence of the 
resolution as a function of total charged track energy in each decay. This is done separately 
for each of the four different b hadron types. Offsets in the boost reconstruction, especially 
for decays with low reconstructed boost, have been corrected for as well. 

The decay length resolution is also parametrized by the sum of two Gaussians. 
For decays with more than 1 track in either the B or D vertex, the resolution <jl is 
estimated from the B vertex fit and IP position measurement errors, combined to yield an 
uncertainty <J meas along the flight direction. The decay length resolution is then obtained 
by appropriately scaling this quantity to determine a 60% core resolution aii — s\ x a meas 
and a 40 % tail resolution <tl2 = $2 x a meas . For correctly tagged B® decays we find 
s\ = 0.92 and s 2 = 2.30. For decays with 1 track in each of the B and D vertices, 
the resolution is extracted from the overall decay length residual distributions in the 
simulation. The difference in decay length resolution between right and wrong charge 
dipole tags motivates treating those separately in the likelihood function, see Eq. (|i~5|). 
(This is similar to differences in resolution between (b — > I) and (b — > c — > I) in the 
lepton+D analysis.) For example, the average decay length resolution for all B° s decays 
with right (wrong) charge dipole tag can be parameterized by the sum of two Gaussians 
of widths cr^i = 76 /iin (112 /im) and o l 2 = 311 fim (450 /im), where the first Gaussian 
represents 60% of the decays. In the case of B® and B° s decays, the right- and wrong- 
tag decay length resolutions are estimated for each of the five decay final states. Decays 
reconstructed within 200 /im of the IP have worse decay length resolution and suffer from 
asymmetric tails (this is most likely due to the addition of a primary track in the secondary 
vertex). These effects have been taken into account in the analysis. Finally, offsets in 
the reconstructed decay length are corrected separately for decays involving one or two 
charm particles, as the effect is small for the former but not negligible for the latter. 

5.4 Oscillation Analysis 



An amplitude fit is performed, as described in Sec. [4.2| and the result is displayed in Fig. |TJ 
Systematic uncertainties are estimated as for the lepton+D analysis except for those 
affecting the final state tag. Here, uncertainties in the charge dipole correct tag probability 
modeling are obtained by varying the scale factor applied to the probability derived from 
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Figure 14: Measured amplitude as a function of Am s in the charge dipole analysis. 



the MC simulation according to So = 0.95±0.025. In addition, the systematic uncertainty 
also includes the effect due to a further reduction of the scale factor by 5% for decays within 
L < 0.5 mm of the IP. Dominant uncertainties are due to the production fraction in 
Z° —>■ bb events, the boost and decay length resolutions, and the overall uncertainty in 
the final state tag purity, see Table B 
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Table 2: Measured values of the oscillation amplitude A with a breakdown of the main 
systematic uncertainties for several Am s values in the charge dipole analysis. 
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-0.004 
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+0.021 
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decay length resolution 


+0.017 
-0.020 


+0.048 
-0.035 


-0.026 
-0.002 


+0.080 
-0.047 


boost resolution 


+0.049 
-0.044 


+0.209 
-0.227 


+0.075 
-0.096 


+0.450 
-0.084 


B® lifetime 


+0.039 
-0.038 


+0.014 
-0.017 


+0.030 
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+0.028 
-0.045 


Am d 


+0.002 
-0.001 


+0.001 
-0.001 


+0.000 
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+0.001 


initial state tag 


+0.001 
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+0.016 
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+0.004 
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final state tag 


+0.069 
-0.099 


+0.036 
-0.236 


+0.027 
-0.110 


+0.074 
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6 Combination of the Analyses 

The D s ±tracks, lepton±D and charge dipole analyses are combined taking into account 
correlated systematic errors. Events shared by two or more analyses are assigned to the 
analysis with the best sensitivity such as to produce statistically independent analyses. 
Figure [3j| shows the measured amplitude as a function of Am s for the combination. As 
noted earlier, the measured values are consistent with A = for the whole range of Am s 
up to 25 ps" 1 and no evidence is found for a preferred value of the mixing frequency. The 
following ranges of B®-B® oscillation frequencies are excluded at 95% C.L.: Am s < 7.6 
ps" 1 and 11.8 < Am s < 14.8 ps" 1 , i.e., the condition A ± 1.645 < 1 is satisfied for 
those values. The combined sensitivity to set a 95% C.L. lower limit is found to be at a 
Am s value of 13.0 ps" 1 . These results are preliminary. 
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Figure 15: Measured amplitude as a function of Am s for the lepton+D, D s +tracks, and 
charge dipole analyses combined. 
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